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Abstract 

We study the block spin transformation for the 2D Ising model 
at the critical temperature T c . We consider the model with the con- 
straint that the total spin in each block is zero. An old argument by 
Cassandro and Gallavotti allows to show that the Gibbs potential for 
the transformed measure is well defined, provided that such model has 
a critical temperature T' c lower than T c . After describing a possible 
rigorous approach to the problem, we present numerical evidence that 
indeed T' c < T c , and a study of the Dobrushin-Shlosman uniqueness 
condition. 
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1 Introduction 



In this note we discuss the block spin transformation for the two dimensional 
Ising model at the critical point, trying to show that it is well defined, and 
it gives rise to a Gibbs measure corresponding to a translationally invariant 
finite norm potential. If one wants to define the renormalization group flow 
in the space of Hamiltonians this is of course an essential point. Here we 
study only the first step of this flow. Similar results have been obtained by 
Kennedy [|T^] for the majority rule renormalization group transformation. 

Our aim would be to use a rigorous approach to the subject but, as it will 
appear clear in the sequel, this is presently a very difficult task. Hence we 
will mainly analyze the problem from a numerical point of view. However, 
as we will explain, this will be partially achieved by combining theoretical 
perfectly rigorous ideas with numerical tools in order to "measure" via a 
computer some interesting theoretical quantities. 

The problem of considering well defined renormalization group transfor- 
mations (RGT) and, in particular, the study of RGT in the framework of 
the modern rigorous approach to statistical mechanics, has attracted the 
attention of several authors. 

Recent and less recent papers have been dedicated to this subject. We 
want to quote in particular, among the "old" papers, the one by Cassandro 
and Gallavotti [[| (which, at our knowledge, is the first one treating explicitly 
the above mentioned problem), the paper || by Griffiths and Pearce and the 
one by Israel fl2[]. For the recent results we quote the monumental paper [|7| 
by van Enter, Fernandez and Sokal, where the problem is discussed in a very 
clear and complete way. This article-book is reach of very interesting dis- 
cussions and examples. It contains a self-contained exposition of the general 
set-up and a very extensive and up-to-date bibliography. Thus we refer the 
interested reader to for a review on the subject and for reference to other 
recent papers. 

It is worthwhile to remark that the majority of the examples considered 
in the above papers, concerning the rigorous approach to RGT, deal with 
the region far (and often very far) from the critical point. The paper 
constitutes, in some sense, an exception and some crucial ideas developed in 
the present paper go back, in fact, to 0. 

In the following section we shall present a general critical discussion about 
the definition of RGT. In section |] we shall give the definitions concern- 
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ing the models considered in this paper and the Monte Carlo procedure 
used to study their equilibrium properties. Section || will be devoted to the 
Dobrushin-Shlosman uniqueness condition and to the exact definition of a 
related numerical quantity, that we have studied to analyze the RGT (see 
next section). Sections || and |^ contain the numerical results and, finally, 
section [7] is devoted to the conclusions. 

2 A Critical Discussion About the RGT 

The main question arising when rigorously discussing the RGT can be ex- 
plained by considering, for example, the Ising model at magnetic field h and 
inverse temperature (3. Let 

v — T(b) fip,h 

be a measure arising from the application of a renormalization group trans- 
formation T(6), defined "on scale b" , to the Gibbs measure //^ of the Ising 
model. Transformations of this kind are always trivially defined in finite vol- 
ume but, of course, we are interested in taking the thermodynamic limit or, 
rather, in defining directly the transformation in the infinite volume situa- 
tion. 

The main "pathology" that can take place, which is the main part of 
the discussions contained in 0, is that v can be non-Gibbsian. This means 
that the conditional probabilities of v can be incompatible with the Gibbs 
prescription corresponding to any absolutely summable potential. This non- 
Gibbsianness is detected via the violation of a necessary condition, namely 
the property of quasi-locality for the conditional probabilities of v. 

The Dobrushin-Lanford-Ruelle (DLR) theory of Gibbs measures is bas- 
ed on the conditional probabilities tt\ for the behavior of the system in a 
finite box A CC Z d , subject to a specific configuration in the complement 
of A (we use the notation A CC Z d to denote a finite subset A of Z d ). For 
simplicity let us only consider Ising-like systems. The configuration space of 
the system, in this case, is Q = { — 1, 1} Z ; we use f2 A = { — 1, 1} A to denote 
the configuration space in A C Z d . According to 0, a probability measure 
whose conditional probabilities {7TA}A C cz d satisfy: 

Hm sup |tta/(wi) -n A f{u 3 )\ = (1) 

A'T Z (Ji,tJ2efi:(an) A /=(aJ2)A' 
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(namely the conditional expectations in A of any cylindrical function / cor- 
responding to different boundary conditions u)i,u)2, coinciding in A' D A, 
tend to coincide as A' tends to Z d ), is called quasilocal. ([I]) can be seen as 
a continuity property in the conditioning, infinite volume, configuration u. 
Kozlov has shown in [14j| that a quasilocal probability measure on Q which 



also satisfies a so called non-nullity condition, i.e. a sort of absence of hard 
core exclusion, is Gibbsian, in the sense that its conditional probabilities 
can be obtained, via the Gibbs prescription, from an absolutely summable 
potential (see 0] for more details). 

It is useful to make at this point some remarks about this notion of 
quasilocality. 

1) Kozlov's theorem (i.e. the fact that nonnullness together with quasilocal- 
ity imply Gibbsianness) is proved in [H| by using an approach which can be 



considered somehow artificial. Starting only from some nice continuity prop- 
erties of the conditional probabilities {vta}, one gets a series, representing the 
interaction of a point x with the rest of the world, which is, a priori, only 
semi- convergent. One can insist to extract the many body potentials from 
that series (for instance via the Moebius inversion formula) pretending that 
they are absolutely summable. This can be achieved only by regrouping the 
terms in some suitable order. One can use, for instance, the lexicographic 
order of the lattice. The resummation will depend in this case on the location 
of x. The potential will be absolutely summable but, in general, not trans- 
lationally invariant. To get translational invariance one needs some stronger 
properties on how weakly the conditional probabilities depend on far apart 
configurations. 

In some situations to compute the renormalized potentials one can use 
much stronger methods, based on convergent cluster expansions; in this way 
a genuine finite norm, translationally invariant, potential is produced in a 
very natural way. Each renormalized coupling constant is expressed via a 
convergent series (see, as an example, the paper JIJ by Cammarota). 

2) The above notion of quasilocality of a measure [i needs a control uni- 
form in the conditioning configuration. To prevent the existence of a Gibbs 
potential it is sufficient that the condition is violated for only one special 
configuration; in this case even the somehow artificial quantity introduced 



in [13] cannot be constructed. However, for any infinite volume reasonable 
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stochastic field /i, a single infinite volume configuration 00 is of zero mea- 
sure; moreover it can even happen (see below) that the single configuration 
inducing non-Gibbsianness, as a consequence of non-quasi-locality, is very 
"non-typical" with respect to /1. It is then natural and physically relevant 
to introduce a weaker notion of quasilocality, for instance by requiring the 
validity of a condition like (D) only for /i-almost all u's. 

A precise definition in this sense has been recently introduced by Fer- 
nandez and Pfister (see ||), but this does not prevent the construction of 
pathological examples. In fact in || the authors show that, for some inter- 
esting examples, a very strong notion of non-quasilocality holds, in the sense 
that (|l|) fails actually for [i-almost all configurations to. 

This happens, for instance, for the example of non-Gibbsianness given 
by Schonmann in p3| . This example consists simply in considering the rela- 



tivization v of the measure /i + for the Ising model in two dimensions to the 
line I — {x = (xi,X2) G Z 2 : x 2 — 0} (isomorphic to the one dimensional 
lattice Z 1 ). Here we are at large inverse temperature (3 and zero magnetic 
field; /i + is one of the two extremal Gibbs measures, the one obtained via a 
thermodynamic limit with + boundary conditions; finally, by relativization 
of /i + to / we simply mean the projection on the <r-algebra generated by the 
spins in I or, simply, the (marginal) distribution of the spins in I deduced by 
ji + by integrating out all the spins in Z 2 \ /. 

3) There are many cases in which a stochastic field /x shows up the pathol- 
ogy of non-Gibbsianness, like in the Schonmann example; however, at the 
same time, the measure fi can be, from many respects, very well behaved. 
For instance the one dimensional measure v of the Schonmann example has 
exponentially decaying truncated correlations. Moreover, as it has been re- 
cently shown by Lorinczi and Vande Velde (see [|15j), this even very strong 
non-Gibbsianness is, in a sense, an unstable property. 

Let us reanalyze the Schonmann example. Suppose that, instead of con- 
sidering the one dimensional sublattice I, one considers a sublattice lb of suf- 
ficiently large spacing b. Namely starting from one integrates out all the 
spins outside the set l b = {x = (xi,x 2 ) G Z 2 : x 2 = 0, X\ = nb,n G Z 1 }, 
obtaining the relativized measure on l^. 2/& can also be seen as obtained 
via a decimation procedure on scale b from v. In general, given a measure 
/i on VI = {-1,1} Z and an integer b > 2, the decimation transformation TJ, 
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acts on fi so that 

v = T b fi 

is simply the relativization of fi to the sublattice Z d of Z d with spacing b, 
that is Z d = {x G Z d : x = by, y G Z d }. Lorinczi and Vande Velde show 
that uj, is Gibbsian in the strong sense, that the renormalized potential can 
be computed via a cluster expansion and it is absolutely convergent. 

Another interesting example given in M concerns the decimation trans- 
formation applied to the unique Gibbs measure /^g^ for the Ising model at 
large f3 and h ^ 0, say h > 0. They show that, V b and for suitable f3 and h, 
the renormalized measure v = Tbfip^h, arising from a decimation transforma- 
tion with spacing b, is not consistent with any quasi-local specification. In 
particular it is not the Gibbs measure for any uniformly convergent interac- 
tion. As it is noticed in |7j] the non existence of the renormalized interaction 
is a consequence of the presence of a first order phase transition for the orig- 
inal model in Z d \ Z d for particular values of (uj x ) x£Z d and suitable values of 
h and (3. One example is the case where uj x = — 1 Vx, h uniform and positive, 
exponentially in (3 near to the value h* (b) which is needed to compensate, in 
Z d \ Z d , the effect of the — l's in Z d and to give rise to a degeneracy in the 
ground state in Z d \ Z£ (the highly nontrivial part in the proof, given in [^], 
of the existence of the pathology consists in showing, via the Pirogov-Sinai 
theory, the persistence of the phase transition at positive temperatures). 

On the other hand from the above analysis it is clear that this pathology 
comes from the fact that, on a too short spatial scale b (with respect to the 
thermodynamic parameters and mainly to the magnetic field h), the system 
is reminiscent of the existence of a phase transition for h = 0. It seems 
reasonable that this pathology could be eliminated provided one uses a RG 
transformation defined on a proper scale depending on the thermodynamic 
parameters. In Martinelli and Olivieri have shown, exactly for the above 
example of Ising model for which in J7| the pathology is found, that, with 
the same values of (3 and h, provided one chooses a sufficiently large spacing 
b' > b, the resulting measure Tyiip^ is Gibbsian in the strong sense and that 
the renormalized potential, which is absolutely summable, can be computed 
via a convergent cluster expansion. In particular, taking b' = b n , with n 
sufficiently large, one shows that, iterating n times the transformation T},, 



one goes back to Gibbsian measures; moreover it has also been shown in [[19 
that T^iip h converges, as n tends to oo, to a trivial fixed point. 
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Let us now describe an example of pathology discussed in |7j which is 
particularly relevant in the context of the present paper. It refers to the 
block averaging transformation (sometimes called Kadanoff transformation). 
Suppose to partition Z 2 into square blocks Bi of side 2 (each block containing 
4 sites). The block averaging transformation TM\ consists, in this case, in the 
following transformation applied to the Gibbs measure \hp^ for the Ising 
model at inverse temperature (3 and magnetic field h; the new measure is 
obtained, starting from the original spin variables a x , by assigning to any 
block Bi an integer value rrii and by computing the probability, with respect 
to the original Gibbs measure fJip t h, of the event X^eB; Cx = m %- 

One obtains, starting from fip^iWx})' 

The original system of a x variables distributed according to fj,^ is called 
object system, whereas the new variables rrii distributed according to v con- 
stitute the image system. 

The pathology in the block averaging transformation for the Ising model 
at large inverse temperature (3 and arbitrary magnetic field h is a consequence 
of the existence of a phase transition for the object system for particular 
values of the image variable m r The authors of |7J show that for the configu- 
ration with rrii = 0, V i, the corresponding object system, a constrained Ising 
model, exhibits a phase transition with long range order. As a consequence 
of this fact they are able to show the violation of the quasi-locality condition. 

Of course, since the local magnetizations m ; in the blocks Bi are fixed 
and all equal to zero, the value of h is totally irrelevant. On the other hand 
if h is very large and, say, positive, the object system without any constraint 
is almost Bernoulli with a high probability to have an individual spin equal 
to +1 and the configuration with m 8 = 0, V i, is expected to be very un- 
likely and, in a sense, irrelevant. Probably the weaker condition of almost 
sure quasi-locality introduced in || is satisfied in that situation. Moreover, 
even though v is not Gibbsian, it could have nice properties and its non- 
Gibbsianness could be unstable with respect to small changes. The situation 
could be similar to the previous mentioned phenomenon discovered by Lor- 
inczi and Vande Velde for the decimation applied to the measure appearing 
in the Schonmann example. 

It is interesting now to discuss in some detail one of the main ideas con- 
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tained in the paper which is, in a sense, at the basis of the present note. 
We will consider the block averaging transformation we have just discussed, 
and apply it to the Ising model Gibbs measure /ig^. The authors of 0] are 
concerned, in particular, with the most interesting example where h = and 
(3 = (3 C , i.e. the system is at the critical point. They show that, at least 
formally, it is possible to compute the renormalized potential and to show 
that it is absolutely summable, provided the constrained Ising model with 
all the rrii = 0, V i, is above its critical temperature. 

To be more precise, let H^ r \{rrii]) be the renormalized Hamiltonian cor- 
responding to the renormalized measure z/({mj}) = TM\fJ>p,h', suppose to ex- 
tract from H^ r \{rrii}) all the many-body potentials $ A({^i} ^a) for any 
finite set A of blocks Bi. The authors of show that all the $^'s can be 
expressed as thermal averages of suitable local observables with respect to 
the Gibbs measure corresponding to an auxiliary intermediate Hamiltonian, 
that we call H^\S). H^ 6 \S) corresponds to the constrained Ising model 
with all the m^s set equal to zero. The new (intermediate) local variables Si, 
defined for any block B iy take values in the finite space, containing six states, 
corresponding to the six spin-cr configurations in Bi such that J2xeBi °~x = 0. 

The starting point of the present note is to try to rigorously show, in a 
strong sense, that the auxiliary model with Hamiltonian H^\S) does not 
undergo a phase transition at (5 — j3 c . A priori there are no reasons, as it 
will be clear from the discussion in the following sections, for the critical 
temperature to decrease after the introduction of additional constraints to a 
spin model. We will show for example that, as a consequence of a remark 
due to Kasteleyn, a particular constrained model obtained from the Ising 
zero field model has exactly the same critical temperature as the original 
Ising model ! 

In particular, to detect the absence of phase transition, we will use the 
idea of exploiting some finite size condition, which goes back to Dobrushin 
and Shlosman (see ||, [5], |6|). The basic point is that if one is able to verify 
a condition involving mixing properties of (finite volume) Gibbs measures in 
a suitable set of finite regions, then one can deduce nice properties (typical 
of the one phase region) for the infinite volume system. That can be done, 
for example, by using a computer. One can show for example uniqueness of 
the infinite volume Gibbs measure, analyticity of the infinite volume ther- 
modynamic and correlation functions and exponential decay of truncated 
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correlations. In j|, ^ [J the authors avoid the use of cluster expansion; in 
[|, §] they use conditions referring to arbitrary shapes. 

In Q the authors introduce a somehow weak condition implying only 
uniqueness of the infinite volume Gibbs state and some decay properties of 
the infinite volume truncated correlations. This condition refers to a region 
V CC Z d and is usually called DSU(V) from Dobrushin-Shlosman unique- 
ness condition (see (|19D , (|20|) below). In f|, [| they treat the so called com- 
pletely analytical interactions, proving, on the basis of a stronger condition, 
much stronger results, in particular uniform analyticity and exponential de- 
cay of truncated correlations for any finite or infinite volume A with constants 
uniform in A. 

In pi] , Olivieri and Picco consider similar finite size conditions but 
only for sufficiently regular regions and get similar results of strong type 
(like Dobrushin-Shlosman complete analyticity) by using a block decimation 
procedure and the theory of the cluster expansion. In a series of papers 



( jig , [T7| , |T8| , |19| ) Martinelli and Olivieri developed a critical analysis of the 
known finite size conditions getting new results both for the equilibrium 
(Gibbs state) and for the non-equilibrium (Glauber dynamics) situation. The 
theory developed in jn| |18[ allows, contrary to the Dobrushin-Shlosman 
analysis, to treat, for quite general lattice systems, almost the whole one- 
phase region (see, in particular, jI7| for more details). 

In a recent paper [fJUJ Martinelli, Olivieri and Schonmann showed that, 



in two dimensions, two finite volume mixing conditions of a priori different 
strength called, respectively, Weak Mixing (WM) and Strong Mixing (SM) 
conditions, are in fact equivalent for sufficiently regular domains A (see pOf 
for more details). In |4[] the authors show that if there exists a a region 
V CC Z d such that their finite size condition DSUiV) is satisfied, then weak 
mixing holds for any finite or infinite A C Z d . Then, combining the results 
in |4[] with the ones in |2(J one gets that in two dimensions, if there exists a 
finite region V CC Z 2 such that DSU (V) is satisfied for the constrained Ising 
system with Hamiltonian H^(S), then, for a large class of regular domains, 
including for instance any cube, the strong mixing condition is satisfied for 
this constrained system. 



From strong mixing, using the results obtained in [J21], [22| , one can easily 
make completely rigorous the above mentioned argument introduced by Cas- 
sandro and Gallavotti, and compute the renormalized potentials as conver- 
gent series via the cluster expansions. In this way after proving DSU(V) the 



8 



Gibbsianness of the renormalized measure would be proven in the strongest 
possible sense. 

3 Definition of the Models and of the Heat 
Bath Dynamics 

In the following we will give precise definitions about some models that we 
are going to discuss: the usual 2D Ising model and some "restricted" models 
obtained from the Ising model by imposing some "extensive" restrictions. 
Suppose to partition Z 2 into 2x2 squared blocks Bi, each containing 

4 sites. Each block Bi can be characterized by the coordinates of its lower 
left-hand site yf, namely £>,, = B Vi , where jji = 2x{, Xi G Z 2 and for 
x = (x«,a; (2) ) e Z 2 : 

B 2x = {z G 7? : 2x U) < z (j) < 2{x® + 1), j = 1,2} 

The formal Hamiltonian associated to the usual Ising model in zero magnetic 
field is given by: 

H isin 9 s _ a x a y , (2) 

<x,y> 

where the sum runs over the pairs of nearest neighbors sites in Z 2 , and 

o x e {-1,+1}. 

In the following we will consider a system enclosed in a finite squared 
region A with various boundary conditions; if not explicitly specified, it will 
be understood that the boundary conditions are periodic. 

In the original Ising model there are, in each block Bi, 16 allowed config- 
urations. Instead of the original variables to describe such configuration 
we can as well use the block variables, say Si G {1, . . . , 16} . In each block 
there will be a self-interaction and the mutual interaction between blocks 
deriving from (0) is again of nearest-neighbor type. The Hamiltonian of the 
Ising model, expressed in terms of Si's block variables, will be denoted by 
H^ W \S) = H Isin9 (a). 

We will consider a modified model in which in each block Bi the sum 
m i — J2xeB z °x is constrained to be zero. Now the block variables S% will 
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assume only 6 different values, corresponding to the following six block con- 
figurations: 



+ + 





- + 








+ - 




+ - 




5 


- + 


5 


+ + 


5 


+ - 


5 


- + 


5 



The corresponding Hamiltonian will be denoted by H^(S). 

It is easy to convince oneself that the model with Hamiltonian H^(S) 
has four periodic ground states (see 0). The first one is given by: 

+ + + + + + + + 

:::::::: (4) 
++++++++ 

The second one is obtained from the first one by interchanging the + with the 
— ; the last two ground states are obtained from the first two by interchanging 
the rows with the columns. Note that the last two block configurations in 
(|3p are quite different from the first four; in fact they are the only ones 
that carry a non-zero internal energy and they are absent from the T = 
ground state structure. We call them turnons, since, as we will see, they 
play an important role by allowing the layered ground states to break and 
mix, destroying long range order. They are 



" + 




' - + " 


- + 


5 


+ - 



(5) 



By further restricting the allowed block configurations, so to forbid the pres- 
ence of turnons, we define a new model whose Hamiltonian will be denoted 
by HW(S). 

In the following we shall denote the three models introduced before by 
l( n \ with n — 4, 6, 16. In section [JTT] we will show that with periodic 
boundary conditions is exactly equivalent to two uncoupled Ising models (in 
a smaller volume); hence its critical temperature is exactly the same as in 

j(16)_ 

We have studied the 3 models 2^ by a Monte Carlo procedure, based 
on a suitable Heat Bath dynamics, whose invariant distribution is the finite 
volume Gibbs measure. This dynamics has been used to compute mean 
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values (with respect to the Gibbs measure) of some relevant observables as 
time averages; the mean value will be denoted by (•) in the following. 

We have built a discrete time Heat Bath dynamics based on locally equi- 
librating the Si block variables. We suppose that A is a cube of even side 
size, so that it can be exactly partitioned into N(A) Bi blocks, ordered in a 
lexicographic way, and we define Vt^ = {1, . . . , n} N ( A \ The dynamics in a 
finite volume A is given by a Markov chain defined below. 

We perform a complete update of all the N(A) block variables by suc- 
cessively updating each one of them. For any i = 1, . . . ,N(A), we choose 
at random the new block configuration in Bi, given the configuration 
{Sk}k=i m A, according to the equilibrium Gibbs measure in Bi with S\b9 
boundary conditions (Bf = A\Bi). The related transition probability in the 
X( n ) model is then given by: 

P (n) K _ q>) = exp[-/^(")(^)] 

This is the most efficient method as far as local updates of block variables 
are concerned (since locally we bring at equilibrium the basic block) and is 
implemented by means of a simple look-up table. It is straightforward to 
build up the dynamics so that one can move from one model to another 
simply changing the number of allowed states (the 6 blocks of (0) are stored 
in the first 6 positions of the tables, with the two turnons in position 5 and 
6). 

In order to characterize the critical point of the system we have computed 
two different quantities. First of all we have considered the specific heat Ca 
as defined from the equilibrium energy fluctuations: 

C A = \A\- l (3\{H 2 ) - (H) 2 ) . (6) 

We have also considered a correlation length £ defined by measuring zero- 
momentum correlation functions. We define the sums over planes 



(*) = E *M , (7) 



s=l 



where x = (r^r 2 ') = (t, s), and the correlations 
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G ^ = Tr~t\ E (™(*o)™(*o + t)) . (8) 

I I to=l 

A t-dependent correlation length (which we will plot for t — 5, where we get 
a fair estimate for its limit as £ — > oo) can be defined by 

For the constrained models, which have a pathological behavior at odd sepa- 
rations, we have found practical to define the correlation length by a distance 
2 ratio 

fe« = (i log ^_)-. do) 

Note that in the T^ 16 ) model coincides with £(2)(0 for t — > oo. 



3.1 4 State Restriction and Full Ising Model: Proof 
of the Equivalence 

Let us denote Bl the family of 2 x 2 blocks partitioning the cube of (even) 
side size L and {S a } ae B L the generic configuration of the 1^ model. We 
have: 

H^\S)= £ w(S a ,Sfj) (11) 

<a,/3> 

where < a,/3 > denotes a couple of nearest neighbor blocks and w(s a , sp) is 
the interaction energy between the two blocks (the self interaction vanishes). 

We want to show that it is possible to associate to each block a an 
invertible map <f> a : S a —> (p a ,T a ) G {— 1,+1} 2 , so that, for any couple 
< a, (3 > and any choice of S a , Sp: 

Vj(S a , Sp) = p a pp + T a Tp (12) 

We first observe that all the four block configurations are of the form 

-t7 2 -01 



(13) 
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Hence there is a simple way to define 4> a on the blocks contained in a 4 x 4 
square, so that (|i~2|) is verified at least for the blocks in the square. This 
definition can be visualized in the following picture: 



-r 2 


-P2 




~P3 


-r 3 


P2 


T 2 




T 3 


P3 


Pi 


T\ 




T 4 


P4 




~Pl _ 






-T 4 _ 



(14) 



Moreover, if L is a multiple of 4, there is a partition of the lattice into 
4x4 squares and it is immediate to check that (|T^) is satisfied for all couples 
< a, (3 >, if 4> a is defined in each 4x4 square as in ([T4]), by periodic extension. 

To complete the proof it is sufficient to observe that (|l|) implies the 
following identity for the partition functions: 

ZxW,M = Z i^,l/2(P) 2 (15) 

The model has four ground states, exactly coinciding with the ones 
ofX (6) (the turnons in this case are absent at T = 0). These four states 
correspond, via the map a , to the four ground states of the two independent 
Ising models (all +1 or all —1 for each one of the two Ising systems). 



4 The Dobrushin-Shlosman Uniqueness 
Condition 

Let us define the variation distance between two probability measures /ii and 
fi 2 on a finite set F[| as: 

Var(/i 1; /i 2 ) = - \Pi(v) ~ ^{y)\ = sup \nt(X) - p 2 {X)\ (16) 

2 yeY XCY 

Given a metric p(-, ■) on Y the Kantorovich - Rubinstein - Ornstein - Vasser- 
stein distance with respect to p between two probability measures /ii, fi 2 on 
Y, that we denote by V p (p\, p 2 ), is defined as 

^ much more general framework can also be considered 
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Vpifi!,^) = jpf p(y>y')»(y,y') ( 17 ) 

where K(ni, /i 2 ) is the set of joint representations of \i\ and /z 2 , namely the set 
of measures on the cartesian product Y xY whose marginals with respect to 
the factors are, respectively, given by p\ and // 2 - This means that, \/B C Y: 

fi(BxY)= Yl riv,l/) = Hi(B) , 

yGBy'eY 

For the particular case 

it is possible to show that V^(/ii,/i 2 ) coincides with the variation distance 
Var(/ii,// 2 )- 

A result by Dobrushin and Shlosman [Q concerning the uniqueness of the 
infinite volume Gibbs measures generalizes previous results by Dobrushin 
based on a "one point condition" on Gibbs conditional distributions (see 

I)- 

Let us consider a spin system on Z d with single spin space S and finite 
range interaction. We generalize in an obvious way the notation introduced in 
section ||. Given a metric p on S, we associate to it a metric pa on Qj^ = S A , 
for any A CC Z d , by defining: 

Pa(S a ,S' a ) = Yp( S ^ S 'x) 

We say that condition DSU P (A, 5) is satisfied if there exist a finite set A CC 
Z d and a 5 > such that the following is true: for any y G d + A (the 
set of points outside A whose spins interact with the spins inside A) there 
is a positive number a y such that, for any couple of boundary conditions 
r, t' G Q c a with t' x = t x , Wx ^ y: 

^PaOa^a) < OiyP(Ty,T y ) , (19) 
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and 



£ a y < 8\A\ . (20) 

y£d+A 

where p\ is the Gibbs measure in A with boundary conditions r outside A. 
We say that DSU(A, 5) is satisfied if (£§ and © hold with p given by flU). 



Theorem 1 (Dobrushin - Shlosman j^j) Let DSU p (A,5) be satisfied for 
some p, A and 5 < 1; then 3 C > 0, 7 > such that condition WM P (A, C, 7) 
holds for every A. 

By WM p (A, C, 7) we mean a particular mixing property of /i^, saying that 
the influence at x G A of a change in the conditioning spins r decays as 

^ e -7dist(a;,9A) 



See |T7]] for a precise definition. 

Theorem [TJ implies, in particular, the uniqueness of infinite volume Gibbs 
measures. Then (|T9|), (|20"D provide an example of finite size condition: one 
supposes that some properties of a finite volume Gibbs measure are true and 
then deduces properties of infinite volume distributions. 

We observe now that 

xGA 

where \i T x is the probability distribution of the spin s x with respect to the 
measure /i A . Hence, if p is given by (|l^), then 

^a(^a^a) > E Var (^>^') • 

It follows that, for given A, DSU (A, S) is certainly not satisfied for any 5 < 1, 
provided that: 

sup £ Var(/4X') > , y yed +A 

t,t' :T x =T^.,\/x^y zgA |0 A I 

This observation will be used in our numerical calculations on model in 
the following way. We try to find a "good" lower bound for the quantity 
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sup sup Var 0£>/£) 

yed+A L T,T':T x =T' x ,\/x^y x&Kl 

by calculating numerically J2xe\ L Var(/^, for a "large" number of choices 
of y, r, t' and by considering the maximum among these numbers, which we 
call A L . Here A L denotes a cube of side size L in Z 2 and each point of A L 
represents a 2 x 2 block; hence corresponds to a cube of side size 2L in 
the original lattice. 

We consider a "reasonable" indication that DSU (A, 5) is satisfied for A 
large enough, if A L \d + A L \/\A L \ is decreasing in L and there exists an L c such 
that 

M L = A L/2 ^ft = A L/2 | < 1 , L>L C . (21) 
|Al/2| L 



5 Numerical Results: the Phase Transition 

We will present here a first set of numerical results, which give numerical 
evidence that the model undergoes a phase transition at a temperature 

The correlation length of the constrained model at the critical point of the 
full Ising model is a finite, reasonably small number, which we estimate. 

This numerical evidence is not meant to constitute a large scale simula- 
tion. We do not try here to estimate critical exponents or to determine with 
high precision the position of critical points (for large scale simulations of the 



2D Ising model see for example |l L] and references therein). The main point 
here is to show in a non-ambiguous manner that the two critical tempera- 
tures are different and that the correlation length of the constrained model 
at the critical point of the Ising model is finite 

Our results have been obtained for cubes containing 400 2 lattice points 
(i.e. the region A^ with L = 200 according to the notation of section ^). In 
the case of model we have performed 210 5 full sweeps (that is full update 
of all lattice sites) of our Heat Bath block algorithm per each value of the 
inverse temperature (3. In the case of and Z^ 16 - 1 models we have used 10 5 
sweeps per each value of (3. We have simulated smaller volumes to check that 
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everything was well compatible with the expected finite size behavior (but 
we will not report in detail about these data). All the runs discussed in this 
Section use periodic boundary conditions. 

Let us start by discussing our simulations for the full 2D Ising model. 
As we said, very large scale simulations exist (pl]|) and the results we are 
presenting here are just meant to set the frame for showing numerically the 
different behavior of the two relevant models. So we have simulated the 
full 2D Ising model in the same conditions that we have used to study the 
constrained X^ model. 

In fig. [I] we show the specific heat of the model. The point closer to crit- 
icality is the one at (3 = .4400 (as is fortunately obvious from the picture). 
Here, as well as in the following, the point size is of the order of magnitude of 
the statistical error, except for the point closer to criticality. The mild (log- 
arithmic) divergence of the specific heat is at an inverse critical temperature 
which we can estimate from our data to be at f3 c = .440 ± .001, in agreement 
with the known exact value. 

In fig. |2| we show the correlation length £ = £(5) (see (|9])), which diverges 
at T c with a critical exponent v — 1. We plot £ for (<L. 

The second set of numerical simulations refer to model X^\ We show 
respectively the specific heat and the correlation length in figs. |3| and f|. 

These two figures strongly suggest that Tf (6) ) < Tf (16) ): the critical 
temperature of the constrained model is smaller than the one of the original 
model. At the critical temperature of the Ising model the constrained model 
does not show a critical behavior. We estimate 

/3f (6) ) = .4775 ± .0025 . (22) 

The correlation length of model X^ is finite at the critical temperature 
of the Ising model; we estimate (see (0) for the definition and note that it 
is defined in units of the original lattice): 

£ (2) =£ (2) (5)~ 11.5 ±0.5 , if/3 = /3f (16) ) (23) 

This is in agreement with the numerical calculations related to the DSU 
condition, that we shall discuss in the following section. In fact we find that 
the inequality (|2l|) of section |] is satisfied for L > L c , with L c ~ 5£. 

As we have already said, we have also simulated the model restricted to 
4 states (by forbidding turnons), which is equivalent to two decoupled Ising 
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models. In figures [5] and ^| we show the specific heat and the correlation 
length; they clearly show criticality at T^ lW ^ = Tj? ia \ Note that the corre- 
lation length is larger than in model X^ l& \ as one expects because the effective 
lattice spacing of the two independent Ising models is half the spacing in the 
original lattice. The ratio between the two correlations lengths should indeed 
be exactly 2, if we could calculate the limit of C(2)(t) as t — ► oo. 

In order to better characterize the nature of the transition for the con- 
strained model 1^ we will present some more data. As discussed in section 
H, at T = the 1^ model has 4 ground states, with broken translational 
symmetry. We define pi, % = 1, . . . , 4 the projection of a given configuration 
over the i — th ground state (we count the number of blocks which fit the 
ground state pattern, and we normalize with the total number of blocks). 
We define: 

P = \(Pl + P2 + P3 + Pi) • (24) 

Let us remind again that the ground states do not include turnons. The 
turnon density is (1 — 4p), and tends to zero for (3 — > oo. In fig. [7| we plot p 
as a function of T for model X^. The density of turnons at (3 = .4407, the 
critical point of the full 2D Ising model, is close to .06, and is determining 
the finite correlation length. 

Spontaneous symmetry breaking is signaled by the non-vanishing, in the 
infinite volume limit, of a defined by 

^E(p t -P) 2 - (25) 
4 i=i 

In the symmetric phase a — > in the infinite volume limit. In the broken 
phase the system aligns in one of the 4 ground states (tunneling, at finite 
volume, among the different layered ground states), and a is non zero. We 
plot cr in fig. The location of the critical point is signaled with high 
precision from the drastic change in a. 
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6 Numerical Results: the Dobrushin-Shlos- 
man Condition 

As discussed at the end of section f| we have used the inequality (PD to 
obtain some reasonable insight about the validity of the Dobrushin-Shlosman 
uniqueness condition, which is too difficult to verify also numerically for 
volumes as large as needed in our case. However, even the task of checking 
numerically the inequality (|2T| ) is indeed a difficult one, with respect to the 
calculations of the previous section. In fact in this case one has to control the 
whole distribution. This is quite difficult, if compared with the simple task 
of computing averages of some observable quantities. The internal energy, 
for example, is peaked close to expectation value; hence, in order to compute 
its expectation value, one just needs to explore a very restricted part of the 
phase space. A computation of M. L (see (pT|) ) demands, on the contrary, as 
we will see, a very large statistics. 

In our simulations we have considered a sequence of cubes of side size L 
and, after selecting L and the number I of Monte Carlo sweeps (we define 
a Monte Carlo sweep as the update of all lattice sites), we have chosen N 
different random boundary condition, by assigning to the boundary sites the 
value + or — with equal probability. We shall call r^., k = 1, . . . , N these 
boundary conditions. 

For each we have selected (randomly) a boundary couple of adjacent 
sites belonging to the same boundary block and we have considered the four 
boundary conditions T 3 kl j = 1, ... ,4, differing from each other only in the 
chosen boundary block, by changing in all possible ways the values of the two 
spins. We have done 4 Monte Carlo runs with the four boundary conditions 
and, during these runs, we have recorded the number of times each block 
variable has visited each of the 6 allowed states, by constructing 



where j ranges over the 4 different values, i ranges over the (L/2) x {L/2) 
blocks and Si ranges over the 6 allowed block configurations. N is normalized 
in such a way that 



(26) 



EiVj(*,5i) = l. 



(27) 
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The quantity of ( |21| ) was then calculated as 



In our simulations we have used L — 8, 16, 32, 64, iV of order 100 and 
values of / ranging from 4000 to 4 10 6 (for the largest value of L). 

In fig. |9| we plot the quantity Ml of (^l|) as a function of the inverse 
square root of the run length, for different lattice sizes and number of iter- 
ations. The straight lines are our best linear fit, which turns out to be the 
right ansatz for the observed behavior. 

The inequality (|2T|) is satisfied only on the 64 2 lattice, and it is violated 
on smaller lattices. Our best fits give 



Ms = 2.10 + , 

M w = 1.83 + ^L, (29) 



M 32 = 1.32 + 
M 6i = 0.58 + 



v 7 ! 
229 

77 

706 

77 



There is indeed a big contribution due to the fact that we are adding a finite 
number of positive random numbers. Only in the limit of large number of 
iterations, I — * oo, this contribution goes to zero. In order to minimize this 
effect we have also run simulations with the same boundary condition r^, for 
each k. The average value of the difference 

{WlN^S^-N'^SJl} 

Z i,Si 

obtained in these conditions has been subtracted from Ml in order to define 
Ml- The contribution we have subtracted from M l goes to zero in the limit 
I — > oo, making Ml a good estimator for the inequality (pT|). We plot Ml 



in fig. [10|. Our best fits for Ml give 
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M 32 



M 64 




(30) 



The constant term coincides, as it should, with the one we get for Ai, with 
very good precision. On the contrary slopes are smaller, indicating that M. 
is, as expected, a better estimator than Ai when using a finite number of 
iterations. 

7 Conclusions 

From the numerical results on specific heat and correlation length, which are 
based on "traditional" numerical methods to detect a (second order) phase 
transition, we can be reasonably sure that, indeed, the critical temperature 
of X^ 6 ' model is strictly less than T^ sm9 . The remark on the isomorphism 
between X^> and the full Ising model Z^ 16 ) shows that one cannot, a pri- 
ori, expect any monotonicity property of critical temperatures in terms of 
imposed constraints. Our numerical results indicate that, by imposing the 
constraint = 0, V Bi, one decreases the critical temperature whereas, by 
enhancing again the constraint via the further elimination of the turnons, 
one increases the critical temperature since it goes back to T^ smg . 

To theoretically analyze this apparently strange behavior it seems useful 
to use some generalized form of Fortuin-Kasteleyn representation of the Ising 
model sufficiently "elastic" to include, in the same set-up of random-cluster 
models, the three models X^ 16 \ X^\ X^ that we have considered (JT^]). 

On the other hand to rigorously implement the Cassandro-Gallavotti pro- 
gram one needs a strong notion of absence of phase transition, namely to 
verify some strong mixing condition SM. 

We say that a Gibbs measure fj,\ in A with r boundary conditions outside 
A satisfies a strong mixing condition SM(A, C, 7) if the influence at x G A 
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of a local change in y G A c of the value of the conditioning spin configuration 
r decays as 



Q e -i\ x -y\ 

Of course SM(A, C, 7) implies WM(A,C, 7) (see section |J); moreover S*M 
(A, C, 7) is interesting when it is valid for a class of volumes A invading Z d 
with C and 7 independent of A. 

As we have said in section ^ in |2(J it is proven that, in two dimensions, 



WM implies SM at least for sufficiently regular regions A. SM would be 
certainly largely sufficient to compute, via the Cassandro-Gallavotti method, 
the renormalized potentials whereas, in general, WM alone will not. 

In section [| we have analyzed a sort of lower bound (involving variation 
distance) for the quantity appearing in the Dobrushin-Shlosman uniqueness 
condition (implying WM and so, since we are in 2D, SM). This is only 
an (almost) necessary condition to be verified in order to satisfy the true 
condition DSU. It only gives an indication in the sense of the possibility to 
verify DSU (which involves Vasserstein distance). 

It is clear that one needs to consider squared regions containing about 30 
x 30 blocks. This rules out, at least with the present time computers, the 
possibility of a computer-assisted proof. Then it would be important to be 
able to find a "Montecarlo" method to "measure" the quantity (Vasserstein 
distance) appearing in DSU. 

It would also be interesting, in general, to find an algorithm, easily im- 
plementable in a computer, to evaluate the Vasserstein distance between two 
Gibbs measures with different boundary conditions. This will be the object 
of further investigations. 
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Figure 1: The specific heat (as computed from energy fluctuations), as a 
function of (3 for the Ising model. 
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Figure 2: The correlation length £ as a function of /3, for the Ising model. 
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Figure 3: As in fig. |l|, but for the 6 block state constrained Ising model, 
J( 6 ). 
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Figure 4: As in fig. but for the 6 block state constrained Ising model, 
J( 6 ). 
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Figure 5: As in fig. [I], but for the 4 block state constrained Ising model, 
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Figure 6: As in fig. but for the 4 block state constrained Ising model, 
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Figure 7: p as a function of (3 for the 6 block state model X^ 6 \ 



32 



0.1 



o 
o 



0.05 



0.42 



0.44 



0.46 



0.48 



0.5 



0.52 



Figure 8: a as a function of (3 for the 6 block state model 
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Figure 9: Ml as a function of the inverse square root of the number of 
sweeps, for different lattice sizes. Straight lines are best linear fits. Empty 
triangles for L — 8, filled triangles for L = 16, crosses for L = 32 and filled 
dots for L = 64. 
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